IGPG-07/3-5, NSF-KITP-07-55 
gl^qc/0703144| 



Dynamical coherent states and physical solutions of 
quantum cosmological bounces 

Martin BojowalcQ 

Institute for Gravitational Physics and Geometry, The Pennsylvania State University, 
104 Davey Lab, University Park, PA 16802, USA 

and 

Kavli Institute for Theoretical Physics, University of California, 
Santa Barbara, CA 93106, USA 



Abstract 

A new model is studied which describes the quantum behavior of transitions 
through an isotropic quantum cosmological bounce in loop quantum cosmology sourced 
by a free and massless scalar field. As an exactly solvable model even at the quantum 
level, it illustrates properties of dynamical coherent states and provides the basis for 
a systematic perturbation theory of loop quantum gravity. The detailed analysis 
is remarkably different from what is known for harmonic oscillator coherent states. 
Results are evaluated with regard to their implications in cosmology, including a 
demonstration that in general quantum fluctuations before and after the bounce are 
unrelated. Thus, even within this solvable model the condition of classicality at late 
times does not imply classicality at early times before the bounce without further 
assumptions. Nevertheless, the quantum state does evolve deterministically through 
the bounce. 



1 Introduction 

An understanding of high curvature regimes of a universe is likely to require a quanti- 
zation of gravity which is non-perturbative and background independent. Background 
independence means that one does not base the theory on pre-existing causal or geomet- 
rical structures because they are to be provided by the quantized gravitational field itself. 
Not surprisingly, many difficulties in this new setting have to be overcome even when the 
aim is "only" to verify that a proposed theory will have the correct semiclassical limit. One 
first has to determine appropriate semiclassical states of the interacting quantum theory of 
gravity. Thus, already in the definition of states in which the classical limit is to be probed 
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one has to face quantum dynamics. Unlike perturbative quantizations on a background, 
no exactly known free vacuum state is available which one could use to determine proper- 
ties of an interacting vacuum state perturbatively. Vacuum or coherent states of general 
interacting theories then have to be constructed anew, which can show properties quite 
different from the well-known (Gaussian) states of free theories or the harmonic oscilla- 
tor. Although it is often assumed, Gaussian states may not capture the right semiclassical 
properties in any given system. They may be assumed as "prepared" initial states, but 
crucial deviations from Gaussian form can occur especially in systems with long evolution 
times, for which cosmology is the example par excellence. 

Quantum gravity is not only an interacting quantum field theory whose interacting 
semiclassical states are to be determined, it also, in general, has no close relation to a free 
quantum field theory as it is often exploited in effective field theories of particle physics. For 
correct predictions it is, first of all, necessary to determine precise states which capture 
semiclassical properties. In this paper, a model, introduced in [JJ, is studied which is 
exactly solvable and includes characteristic effects of loop quantum gravity, one candidate 
for a background independent quantization [2J [31 H]. The model itself is based on loop 
quantum cosmology [5]. With new techniques [HI UJ, coherent state properties can be 
determined explicitly. In this sense, the model is analogous to the harmonic oscillator in 
quantum mechanics and it has indeed the same solvability properties as explained in more 
detail below. This will allow us to perform a complete dynamical coherent state analysis, 
demonstrating how properties can differ considerably for distinct systems even when one 
considers only solvable models. The model we study is not only illuminating in this regard, 
but it also is of direct physical interest since it describes non-singular cosmological bounce 
models. 

Bouncing solutions of cosmological models have recently received much attention. Al- 
though they are generally very special, they can indicate how transitions through the 
classical big bang singularity may be possible. Many different examples exist by now, 
which have most systematically been developed in loop quantum cosmology. Most argu- 
ments are based on "effective" equations which import some quantum effects into classical 
equations, and which sometimes allow exact analytical solutions (such as in [El El E]) or 
can at least be studied numerically; see e.g. PU EH E21 E31 El]- Also in this context the 
above question of what a semiclassical state of an interacting quantum theory looks like 
is relevant, although it is often overlooked. It enters in the derivation or justification of 
those effective equations which are supposed to capture properties of semiclassical states. 
If the correct type of semiclassical states is not known, one cannot be sure to have included 
all relevant corrections to the classical equations in the right way. As a byproduct, our 
solvable model presents the first case of a complete set of effective equations in quantum 
cosmology. 

Our model, used to illustrate semiclassical state issues, is solvable exactly at the quan- 
tum level pQ. This is much stronger than having analytical solutions of effective equations 
since full states, including not only expectation values but also fluctuations and higher 
moments, are under full control. In this sense, the system is comparable to a harmonic 
oscillator. A complete analysis becomes possible, including e.g. the evolution of coherent 
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states and their long-term dynamical properties. This is especially relevant in the light of 
recent numerical analyses of related models where initial Gaussian states without squeezing 
were evolved [15] . We will see that the coherent state structure of the model is much richer 
than that of unsqueezed Gaussian states, with squeezing influencing the general behavior 
significantly. This is an instructive example for the importance of a dedicated coherent 
state analysis, rather than taking over harmonic oscillator properties to a new model. 

The model is paradigmatic for background independent quantum gravity obtained from 
a loop quantization where the usual free field theory basis is not available. Since the 
bounce model is exactly solvable, it can provide a perturbative basis for quantum gravity 
including all possible interactions and degrees of freedom. Thus, the form of coherent 
states determined here is relevant not only for the model itself but for quantum gravity in 
general. At this stage, perturbative inhomogeneities are not included explicitly and thus 
the question of how they evolve through a bounce is not addressed in this paper. We rather 
show and emphasize that even the unperturbed isotropic situation poses several important 
questions for how quantum fluctuations of the isotropic mode evolve through a bounce. We 
follow a general method, summarized in Sec. [2j The solvable models relevant for cosmology 
are introduced and analyzed in Sec. [31 and discussed more broadly in Sec. HI 

2 The method 

In what follows, we will not use a fixed, or any, representation of our quantum system on 
a specific Hilbert space. Rather, we take an algebraic viewpoint and treat the algebra of 
basic operators, such as [q,p] = ih in quantum mechanics, together with the Hamiltonian 
H as primary. The quantities we are most interested in are expectation values (q)^ = (q), 
(P)il> = (P) i n a given state ip, which we often drop as a label if no confusion can arise, and 
fluctuations and correlations Aq = a/ (q 2 ) — (q) 2 , Ap = \/ (p 2 ) — (p) 2 and C qp = ~ (qp + 
pq) — (q){p)- Higher moments, involving higher than quadratic powers of basic operators, 
could be included by the same means although we will not need to do so here. Nevertheless, 
it is important to note that, would we determine all the moments, we could reconstruct 
the state ip provided that the moments satisfy appropriate conditions. The most basic 
such condition is Heisenberg's uncertainty relation (Aq) 2 (Ap) 2 — C 2 p > h 2 /4: (also called 
Schrodinger-Robertson uncertainty relation in this form) which will be used frequently 
below. As in [6], we call all fluctuations, correlations and higher moments quantum variables 
since, unlike expectation values, they describe typical quantum properties. 

Determining the evolution of moments is thus sufficient to find properties of states. 
One can side-step the explicit construction of states in a representation by deriving and 
solving equations of motion for moments directly, such as d(p)/dt = {\p,H])/ih and 
d(Ap) 2 /dt = {\p 2 ,H])/ih — 2(p)d(p) /dt. In general, this set of equations is highly cou- 
pled because, unless [p, H] is linear in basic operators, ([p, H]) is a function of expectation 
values and quantum variables. The quantum variables, in turn, will satisfy equations of 
motion involving moments of higher degree. (See [6j [16] for examples.) This coupling de- 
scribes the back-reaction of spreading and deformations on the peak trajectory of a state, 
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which is a crucial quantum effect. It is, for instance, the reason for the usual non-locality 
in time of effective actions. 

This is the place where solvability properties of a model become important. If H is 
quadratic, for instance, [p, H] and [q, H] will be linear in q and p and (q) and (p) will not 
couple to quantum variables. As is well-known for the harmonic oscillator, the spreading of 
states then does not influence the peak motion at all and no non-trivial quantum corrections 
arise in effective equations. More generally, this happens whenever basic operators taken 
together with the Hamiltonian form a linear commutator algebra. Our solvable bounce 
model is precisely of this type. 

It is then feasible to solve equations of motion for (q), (p) and any desired quantum 
variables directly without taking the detour of first computing a state in a chosen rep- 
resentation. Many representation dependent difficulties can be avoided, such as explicit 
formulas for inner products and normalizations^] Instead, properties of the Hilbert space 
structure, such as the self-adjointness of operators, can be implemented straightforwardly 
through reality conditions for the solutions of (q), {p) and quantum variables. 

3 A solvable bounce model and its properties 

A free isotropic scalar field <fi couples to gravity through the Friedmann equation 

We use canonical variables, explained in more detail in [T7j, whose relation to the scale 
factor a, the scalar and their derivatives in proper time r is c = da/dr (extrinsic curva- 
ture), p 3 / 2 = a 3 (volume) and p^ = p 3 / 2 d0/dr. Geometrically, p is the isotropic component 
of a densitized triad and can be positive or negative for the two triad representations. In 
what follows, we assume positive p without loss of generality and drop absolute values. 
Moreover, for later convenience we rescale <p by V% and drop factors of 87rG/3 such that 
{c, p} = 1. Solving (CQ) for p^, which is a constant of motion, yields p v oc ±cp, allowing 
four possible choices for the signs: c can be positive (expanding universe) or negative (con- 
tracting universe), and for each case p^ can take any sign (such that <fi runs opposite or 
along coordinate time r). We can interpret H = cp as the Hamiltonian which generates the 
flow in the variable 0, playing the role of internal time. The deparameterized Hamiltonian 
constraint then reads p^ + H = 0. 

This Hamiltonian is quadratic, although not of the harmonic oscillator form, and thus 
solvable as explained above. Unlike for the harmonic oscillator, the Hamiltonian does 
not have a definite sign. One can easily understand the behavior of solutions, and of 

1 Still, the existence proof of a representation poses an important and sometimes challenging question. 
An explicit representation may also be necessary for foundational issues such as the measurement problem. 
Nonetheless, dynamics can be analyzed without all details of the representation, and is thus useful to keep 
separate from foundational issues which are not fully understood even for quantum mechanics, let alone 
quantum gravity. 
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Figure 1: Illustration of positive and negative energy solutions in an upside-down har- 
monic oscillator potential V(q) = ~\q 2 - The sign of the corresponding triad variable p is 
indicated. 



the energy spectrum of the quantum theory, by performing a canonical transformation to 
new canonical variables (ir, q) in which c = ~^( 7r + ?)> P — ~ q)- The Hamiltonian 
then becomes the upside-down harmonic oscillator H = ^tt 2 — ^q 2 which obviously allows 
positive as well as negative energy solutions. Classical solutions can easily be determined 
as 7r(0) = A cosh + Ssinh0, g(0) = Z?cosh0 + Asinh0. In terms of the integration 
constants A and B, the Hamiltonian is H = \{A 2 — B 2 ). Corresponding solutions in the 
original variables are c(0) = -^(A + _B)e^, p(4>) = -^(A — B)e~^ '. 

Since we assume positive p, solutions as functions of q can only be incoming from the 
left of the upside-down potential, where ir is positive and q negative. This assumption 
implies A — B > 0. The sign of H then depends on whether we describe an expanding 
or contracting universe, c > implying H = | ( A + B) ( A — B) > while c < implies 
H < 0. In the first case, p^ < and (p runs opposite to proper time, while it runs along 
proper time for a collapsing universe. The opposite case would be realized had we chosen 
the opposite sign for H. The behavior of solutions is illustrated in Fig. [TJ 

This Hamiltonian is directly quantized, following the rules of quantum mechanics, in the 
Wheeler-DeWitt approach [TSldH]. It is thus helpful to recall the properties of the resulting 
quantum system as derived in [TJ using the method summarized in Sec. |5J The quadratic 
Hamiltonian implies that the system is solvable, which can be exploited to determine 
explicit solutions (c)(0) = C\t^ and (p)((fi) = c 2 e~^ for expectation values of the basic 
operators c and p, in full agreement with the classical solutions. Although the functional 
expressions for expectation values agree with classical solutions, a difference to the classical 
case is that states may be superpositions of positive and negative eigenfunctions of H even 
though the expectation value of p is either contracting or expanding. We will later discuss 
the effect of possible admixtures of negative to positive energy states in solutions. 
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Moreover, fluctuations can be solved explicitly by the same methods. They are (Ac)(0) 2 = 
c 3 e 2< ^, (Ap)((p) 2 = c 5 e~ 2 ^ and C qp {<p) = c 4 for the correlation. While fluctuations are not 
constant, they satisfy (Ac) /(c) = const and (Ap)/(p) = const. Uncertainty relations re- 
quire C3C5 — C4 > \h 2 . Typically, Ac and Ap are thus of the order y/h, although the precise 
value requires more specific information about the state. For dimensional reasons, this 
must involve the only available scale H such as in the form Ap = \JfiB. and Ac = y^h/H. 
In contrast to the harmonic oscillator where the ground state provides specific values for 
fluctuations, no distinguished state is available for the Hamiltonian encountered here. Later 
on we will see that the loop quantization does provide more information which one can 
use to estimate fluctuations. Still, the information will not be as complete as it is for the 
harmonic oscillator ground state. 

The key feature of a loop quantization is that its representation does not provide a 
c-operator but only operators for almost periodic functions of c [20] (i.e. countable super- 
positions of exp(i/xc) with // G R). This is sufficient for constructing a Hamiltonian operator 
which is well-defined and agrees with the classical one, cp, in low curvature regimes where 
c< 1. A loop quantization can thus be understood as implying that c in the classical 
Hamiltonian is replaced by sin c when it is quantized. Instead of sin c one could, of course, 
choose any other almost periodic function which reduces to c when c< 1. The freedom can 
be restricted by relating the model to quantizations of the Hamiltonian constraint as it can 
be defined in the full theory J2TJ. This distinguishes the use of sine as it follows from the 
original derivations in [17], although other choices remain [22] • In any case, the resulting 
Hamiltonian is no longer quadratic in the variables (c,p) and thus appears quite compli- 
cated. However, if we introduce the operators p and J := pexp(ic), the linear Hamiltonian 
operator H = — ~z( J — .P) reduces to the classical Hamiltonian when curvature c is small, 
and it depends on c only through exp(ic) which is required by loop quantum cosmology. 
It can be viewed as a quantization of p sin c in a specific factor ordering. 

A linear Hamiltonian usually simplifies the dynamics very much, but only if the basic 
variables in which it is linear form a closed algebra (see [6] for a discussion in the general 
context of effective equations). Since we transformed from canonical variables (c,p) to 
non-canonical ones (p, J), the system is not guaranteed to simplify even with a linear H. 
Fortunately, as one can easily check the variables do satisfy the closed sl(2,R) algebra 

[p, J] = hJ , [p, J f ] = -h.P , [J, .P] = -2hp - h 2 . (2) 

which includes the Hamiltonian as a linear combination of J and This is the reason 
for the exact solvability of the system which we will use below. Representations of sl(2, M) 

2 This algebra has repeatedly occurred in canonical quantum gravity, for instance to quantize black hole 
phase spaces [23, 24, 25 . The analog of p in this case is the mass of the black hole, which one requires to 
be positive, selecting the discrete positive series of sl(2, R)-representations. Note the resemblence of the 
basic variables p and J to the variables used in affine quantum gravity [26l [27] . In that case, as with J 
here, one multiplies the metric with momenta to ensure the existence of representations of the resulting 
affine commutation relations where the metric is positive definite. Although we use triad variables and do 
not need to require positivity for them, there are representations where p is indeed positive definite. Also 
in quantum optics, and even for the harmonic oscillator itself, the sl(2,R) algebra plays a role [28, 29, 30 . 
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do not allow one to have a purely positive spectrum for i(J—ji). Thus, any representation 
space contains positive and negative energy solutions and we will have to dicsuss appro- 
priate restrictions on states to rule out superpositions of negative energy contributions to 
positive energy states if superpositions of expanding and collapsing universe branches are 
not to be allowed. 

Before analyzing equations of motion we note that the system can be generalized by 
performing a canonical transformation such as ir := p k c, v := p 1 ~ k /(l — k) and using the 
new canonical variables (tt, v) instead of (c,p) in the definition of p and J. Properties of the 
system to be discussed below are not affected by this re-interpretation of the variables. Such 
transformations can be motivated by taking into account features of nearly isotropic but 
inhomogeneous lattice states of loop quantum gravity and possible dynamical refinements 
of the underlying lattice [HI [32]. The special value k = then corresponds to a fixed lattice 
as realized in [HI [20], while k = —1/2, introduced independently in [32], corresponds to 
a refinement such that the number of lattice vertices increases linearly as a function of 
volume. 

3.1 Equations of motion 

As discussed before, the system can much more easily be understood if we do not first 
solve for wave functions, subject to itvip = Hip, and then compute expectation values 
and fluctuations from solutions. Instead, using a more algebraic point of view we derive 
equations of motion directly for expectation values, fluctuations and higher moments and 
solve them [6l [16] . 

For expectation values, now simply denoted as p := (ft), J := (J) and J := (J') in an 
arbitrary normalized state, equations of motion follow immediately by taking expectation 
values of Heisenberg equations of motion, or by using the Schrodinger equation for the 
state appearing in the expectation value, 

P=±{\p,H]) = -~(J + J) , j=l([j 7 H]) = -p-~h=J. (3) 

That these equations form a closed system is a consequence of the linear nature of the 
variables and Hamiltonian. In general, the evolution of expectation values would also 
depend on all higher moments of the state: During evolution the state spreads and deforms, 
which then back-reacts on the peak position of a wave packet. This back-reaction is the 
dynamical essence of a quantum system, captured in effective equations. 

In fact, a state is characterized not just by its expectation values of basic operators 
but also by the infinitely many quantum variables which specify fluctuations and higher 
moments of a state. Fluctuations (and correlations) can be defined, as usually, by 

qvv := (p 2 )_p 2 (4) 
G JJ := (J 2 ) - J 2 , G JJ := ( J t2 > - J 2 (5) 
G pJ := l -{pJ + Jp)-pJ , G pJ := l -{p.P + J ] p) ~PJ (6) 
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G JJ := \(JP + J ] J) - \J\ 2 - (7) 

Since we use partially complex variables there are initially more than three independent 
fluctuations. However, reality conditions to be imposed later at the quantum level will lead 
to additional relations and reduce the number of independent degrees of freedom to the 
correct value. Higher moments are defined analogously using totally symmetric orderings 
in expressions where both p and J are involved. 

Fluctuations are not expectation values of a single operator and their equations of 
motion do not follow directly as before. But they can easily be derived using linearity and 
the Leibniz rule. We have, e.g., 

qpp = _L ^ h}) - 2pp = -UpJ + Jp + pP + Jtp) + p(J + J) (8) 
in 1 2 

= -G pJ -G pJ (9) 

and similarly 

G JJ = -2G pJ , G JJ =-2G pJ (10) 

1 1 - . - i - 1 

qpJ - _qJJ _qJJ qpp G p ^ —G^^ —G^^ qpp (11) 

2 2 ' 2 2 1 1 

QJJ = _QPJ_QPJ_ ( 12 ) 

Higher moments also are subject to equations of motion which follow analogously. 

For our solvable system, all these equations of motion are linear in the quantum vari- 
ables and only finitely many ones are coupled to each other. They can thus be solved 
straightforwardly, such as 

p{<j>) = \{Ae^ + Be^) - \h, (13) 
J(0) = KAe- 4, - Be*) + iH (14) 

(using —\i{J — J) = H := (H)) for the expectation values. To decouple the six equations 
for fluctuations, we first note that G JJ = G pp and \(G JJ + G JJ ) = G pp implies 

QJJ_QPP = Cl ( 15 ) 

-G JJ + -G JJ - G pp = c 2 (16) 

with constants c\ and Ci- Moreover, only two sets of two coupled equations, one for G pp 
and G pJ + G pJ , 

qpp - G P ^ G p ^ 

^-(G pJ + G pJ ) = - Cl -c 2 -AG pp 
at 
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and one for G pJ - G pJ and G JJ - G 



d 

dt 
d 

d*' 



(G pJ - G p - 



jj 




2(G pJ - G pJ ) 



i(G JJ -G JJ ) 



remain. They yield 



G pp (<P) = \{c,e 



2 " + c^)-\{ Cl + c 2 ) 

2<t> + c 4 e 2 *) + i(3c 2 - ci) - t(c 5 e* - c 6 e"*) 

24> + c 4 e 2 *) + i(3c 2 - ci) + 2(c 5 e^ - c 6 e^) 

2 * - c 4 e 2 *) + l -(c 5 e* + c e e-*) 

*+ - c^e 2 *) - fae* + cee'*) 

2 * + c 1 e 2 *) + i(3 Cl - Ca ) 



(17) 



(18) 



C JJ = ^(c 3 e 

G pJ (0) = \{c,e 

G pJ (<P) = \{c,e 

G J \<P) = \{c,e 



(19) 



(20) 



(21) 



(22) 



with further constants of integration. 

3.2 Physical inner product properties 

Although we are not dealing explicitly with states, properties of the Hilbert space must be 
reflected in the structure of our quantum variables. The equations of motion, for instance, 
are only valid if we understand the quantum variables to be defined with normalized 
states since (|) = 1 has been used, for instance when taking an expectation value of 
[ J, jt] = —2hp — h 2 to compute J in The system is, initially, defined through the 
Friedmann equation as a constraint on the space where both the metric and scalar field 
degrees of freedom are quantized. After quantizing this equation, it is re-interpreted as 
describing relational motion in 0. Loop quantum cosmology provides the kinematical 
inner product on the original Hilbert space where the constraint operator is defined, but 
not immediately one on the solution space. Such a physical inner product can be difficult 
to determine explicitly in a representation of states. It can be derived, for instance, by 
requiring that basic operators are self-adjoint if they correspond to real classical variables. 
This then directly implies that all quantum variables defined from self-adjoint operators 
must be real, which in our procedure is the analog of using the physical inner product in 
deriving expectation values through states. Since such reality conditions can be imposed 
directly for quantum variables, implementing physical inner product properties at this 
level can be much more straightforward than at the level of states. One reason is the 
representation independence of the formalism which allows one to avoid looking for a 
representation of states in which a computation of the physical inner product may be 
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feasible. This is especially useful for quantum gravity where the general physical inner 
product problem is one of the major issues. 

3.2.1 Reality conditions 

In our case, we use one complex classical variable J = pe lc , and thus cannot refer to a self- 
adjoint quantization of c directly since no such operator exists at all in a loop quantization. 
Reality conditions implementing the physical inner product must be formulated in a more 
complicated way: In addition to the simple adjointness relation p^ = p quantizing the real 
variable p, we have a non-linear relation 

JJ ] = p 2 (23) 

which follows from the fact that e tc must be quantized to a unitary operator for c to be 
real. (However, JU ^ p 2 in the ordering chosen for the definition of J = pe lc .) 

Taking expectation values of this equation, we obtain a relation between quantum 
variables and expectation values: using the commutation relations ([2]) in 

G JJ = \(J.P + Jtj) - \J\ 2 = (Jjt) +Hp+ \h 2 - | J| 2 (24) 

we have 

(J.P) = g jJ -h P -^h 2 + \J\ 2 

= (p 2 ) = G pp + p 2 (25) 

and thus 

\J\ 2 ~{P + W = G pp - G JJ + hi 2 = hi 2 - Cl . (26) 

This condition mixes expectation values and quantum variables, but depends on fluctu- 
ations only through the constant c\. An initial state thus determines how the reality 
conditions between expectation values are realized. The relation (1261) is then preserved in 
time since G JJ — G pp = C\ is constant as derived in ffl~5|) . If fluctuations are small, as in 
semiclassical states, we have \ J\ 2 = p 2 + 0(h) which, as it should, is the classical relation 
satisfied up to quantum corrections. 

3.2.2 Sign of the energy 

Although not directly related to the physical inner product, we include in this section a 
discussion of the requirement of a definite sign for the Hamiltonian in superposed states. 
One may or may not wish to allow superpositions of expanding and contracting universe 
branches in quantum cosmology, but arguments using the physical inner product in analogy 
to the Klein-Gordon equation suggest that only energy eigenstates of a definite sign should 
be allowed in superpositions [15]. For a linear system, we did not explicitly take an absolute 
value of the Hamiltonian operator, which implies that in general we are not guaranteed that 
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only positive energy solutions enter states corresponding to our solutions. In a language 
more suitable to quantum cosmology, "positive energy" means that at any fixed time </> 
we should not allow superpositions of expanding and contracting branches of a universe 
(while at different times the universe certainly does not need to be always expanding or 
always contracting). To rule out significant contributions from a superposed branch we 
have to pose further conditions for our quantum variables ensuring that they arise from 
expectation values in states which are superpositions of only positive energy eigenstates 
(or only negative energy eigenstates). In general, expressing the positivity of operators 
through expectation values can be complicated. 

But for our purposes it is, fortunately, possible to proceed without technical com- 
plications. We will be interested in states which at some point (e.g. at late times) are 
semiclassical. This restricts the values that fluctuations can take compared to the magni- 
tude of expectation values. It refers in particular to p^ as one of the matter variables. We 
require that its fluctuation is small compared to its expectation value which, through the 
dynamical equation, implies the relation 

qHH ^} _ (#)2 (#)2 _ (27) 

If this is realized and H = (H) > 0, a state in the ^-representation, i.e. written as a 
superposition of if-eigenstates, is sharply peaked at a large positive value of H. Thus, 
there are no significant contributions from negative energy states. (We will return to this 
issue in Sec. 13.4.31 ) Since H and G HH are constant during evolution, imposing fl27|) at one 
initial time ensures that it is satisfied at all times. 

We can express this condition in terms of the integration constants derived before. 
From 

g hh = ^ 2 )-H 2 (28) 
= -i(J 2 - J J 1 " - J 1 " J + J t2 ) + ^(J 2 - 2\J\ 2 + J 2 ) (29) 
= -\(G JJ + G JJ ) + \g jJ = l -{ Cl -c 2 ) (30) 

we see that C\ — c 2 must be small compared to H 2 . (But it cannot be zero since H 
has continuous spectrum.) This is the primary condition we have to impose not only for 
semiclassicality but also to ensure that a state is dominated by contributions of definite 
energy sign. Later, we will add further semiclassicality conditions to restrict also the 
fluctuations of other variables such as p. 
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3.3 Bouncing solutions 

Our general solution (JT3l) for p allows bouncings! solutions for AB > as well as "singular" 
solutions for AB < which reach p = in finite time (f>. (Although isotropic loop quantum 
cosmology is non-singular for any solution [33J H3], additional correction terms become 
manifest at small volume which are not included here in the solvable model. The model 
itself thus breaks down before p = is reached. The singularity in our equations only 
indicates that a deep quantum geometry regime is reached, just as one commonly expects 
the general singularity problem to be resolved. Nevertheless, we keep solutions with AB < 
for now since they will be ruled out even within our model shortly.) The internal time 
variable has just been chosen for convenience of the mathematical description, rather than 
referring to physical observers. For a solution reaching p = to be considered singular one 
must also verify that proper time remains finite. We thus need to interpret our relational 
solution {p and J as functions of <p) as a space-time geometry subject to modified dynamics 
as it arises from the loop quantization. 

We do not have any manifold picture, except for the homogeneous spatial manifold we 
started with to reduce the classical system. What is missing is a manifold for the time 
extension, which is indispensable if we want to compute a proper time interval. A time 
direction and coordinate can be introduced by reverting back to the Friedmann formulation 
of constrained dynamics. We interpret the effective Hamiltonian density p~ 3 ^ 2 (H) 2 = 
£T 3 / 2 (ImJ) 2 = v /psin 2 c together with the matter contribution ^p~ 3 ^ 2 pi as an effective 
constrainto 

C = -^sin 2 c + ^-p~ 3/2 pl 

which generates coordinate evolution (in if) in a gauge specified by the lapse function N, 
dp/df] = {p, NC}. For proper time, rj = r, we simply have N = 1 and thus 

dp/dr = {p, C} = ^sm(2c) . 

From the equation 

sm(2c(0)j = - — - — = — — (31) 

and our solutions (|T"3l and (TH1) we can then compute 0(r) by integrating 
d0 - 2sin(c(0))cos(c(0)) _ -2V2H 

dr VP[(P> dp/dcj) (Ae~* + Be* - hf' 2 1 } 



3 The bounce occurs at p « H which can be much larger than H although it is quantum effects which 
cause it. Deviations from classical behavior occur in quantum geometry due to large curvature which, for 
a large matter content, can be realized even when the spatial volume is still large compared to the Planck 
length. 

4 We emphasize that this way to derive an effective constraint is correct only due to the linearity in J 
of the Hamiltonian H. Non-linear terms in the Hamiltonian, as they arise from any deviation from the 
solvable model such as by including a matter potential, imply additional correction terms depending on 
quantum variables. They would have to be studied carefully for the derivation of effective equations of 
non-solvable systems; see e.g. [6l [3 [16] for the low energy effective action as an example. 
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(using J/p = cos c + 2 sin c) . We can always assume that either A = B > (for a bouncing 
solution) or A = — B > (for a non-bouncing one) since we only need to shift the origin 
of (j) if \A\ ^ \B\. This leaves us with two cases, 

43/2 r <p 

r (0) = -iL_ / CO sh 3/2 {z)dz (33) 

for A = B and 

A 3/2 /•</> 

r(0) = — / smh 3/2 (z)dz (34) 

for A = — B. The integrals can be determined in terms of elliptic functions, but we are 
only interested in the fact that r(0) is finite at any finite value of which can be seen 
directly from the integrals. Thus, proper time remains finite when p = is reached. 

Singular solutions could thus be possible. But not all these solutions satisfy the reality 
condition (1261) which still has to be imposed. From its general form we obtain 

I J| 2 - (p + \hf = -ab + H 2 = ^h 2 - Cl 

and thus 

AB = H 2 + ci - hi 2 . (35) 

For macroscopic values of H and small (or positive) C\ from fluctuations, we only have 
bouncing solutions with AB > 0. Singular solutions can only be obtained for large and 
negative C\ which is never realized for states which are semiclassical at one time. Note, 
however, that c\ can be large even if our condition G HH H 2 which is necessary for 
solutions to respect positivity is realized since the latter condition only constrains c\ — c^. 
Our discussion thus shows that it is crucial to know and use the reality conditions, or 
ultimately the physical inner product, to draw conclusions about bouncing solutions versus 
non-bouncing ones. 



3.4 Uncertainty relations 

Fluctuations cannot take arbitrary values but are restricted by uncertainty relations. For 
each pair of self-adjoint basic operators we have one uncertainty relation, which in our case 
implies three relations since J + .P and i(J — jr) are independent in addition to p. For 
each pair (A, B) of self-adjoint operators, as usually, we start from the Schwarz inequality 

applied tcfl := AA\if>) and |V> 2 ) := AB\ip) with AA :— A — (A): 

((AA) 2 } ((AB) 2 } > I (A~AAB}\ 2 . (36) 

5 AA is not a linear operator due to the dependence on the state in (A), but \ipi) is well-defined as a 
state obtained from 
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Writing 



and 



AAAB = -(AAAB + ABAA) + i—[AA, AB] 
2 2i 



-(AAAB + ABAA) = -{AB + BA) - AB 
2 2 



we have 



[AA,AB] = [A, B] 



\{AAAB)\ 2 = \{{AB + BAf - 2{A){B)f + ±{-i[A, B]) 2 



where we used self-adjointness of the operators to compute the absolute square of the 
complex number {AAAB). In terms of quantum variables, we thus have the general form 

g aa g bb _ ^ G ABy > ^_i[A,B]) 2 (37) 

of uncertainty relations whenever A and B are self-adjoint. 

Specifically, we have three pairs (p, J + Jt), (p,i(J — ir)) and (J + 3\i{3 — J" 1 ")) of 
different self-adjoint basic operators. We obtain three inequalities involving G pp and the 
fluctuations 

g j + j,j+j := + _py) - {j + jf 

= G JJ + 2G JJ + G JJ = 4G PP + 2( Cl + c 2 ) , (38) 
Gi (j-J),i(j-J) ._ j]y) + {j- jf 

= AG HH = 2(ci — C2) , (39) 
G p ' J+J := \{p{J + J f ) + {J + P)p) - p{J + J) 

= G pJ + G pJ , (40) 

a**'-* := \w - i f ) + (J - i f )p) - < P (J - J) 

= i{G pJ -G pJ ), (41) 
i 
2 



g j+JMJ-J) := 1((J + jt)(j _ jt) + (J - jt)(j + jt)} _ j(j + J)(j _ J) 



= i(G JJ - G JJ ) . (42) 
Using the explicit solutions, they are 

G pp G j + j,j+j _ (Qp,J+jy = Ac^- l -{ Cl + c 2 f (43) 



> Irll 1 



G PP G i(j-J),i(J-J) _ ( G P,i{J-j)f = ( Cl _ C2 )( C3e -20 + C4e 2^ + I( c 2_ c 2) (44) 
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-cje 2 * - 2c 5 c 6 - cle- 2 * 

> + Jf = h 2 (A 2 e~ 24, - 2AB + B 2 e 2 *) 

G J + J,J + J G i { J-J ) , ii J-J ) _ ( G J + J,i { J-J)y = 4(ci _ C2)(c3e -2^ + C4e 20 ) _ 2(c 2_ c 2 ) (45) 

-4c'e 2 * + 8c 5 c 6 - Acle- 2 * 

> h 2 (2p + hf = h 2 (A 2 e- 2,l> + 2AB + B 2 e 2 *) . 

These uncertainty relations are the equations which determine properties of semiclassi- 
cal, near coherent states. We will later discuss these relations, in particular their saturation, 
in more detail and give a complete analysis of coherent states of this system. Before doing 
so we can already note here that there are quite unfamiliar properties compared to what 
one knows from harmonic oscillator coherent states. The first relation, (T4"5|) shows that 
there is a type of uncertainty relation between the constants of integration c 3 and C4, one 
of which determines the p-fluctuation before and one after the bounce. Thus, if the un- 
certainties are very small at very late times, say, they must have been very large at early 
times. This relation also indicates that equally distributed fluctuations are typically of the 
size ~ VhH and thus (Ap) /p = VG^/p ~ \fhfU. 



3.4.1 Saturation 

Of particular interest is the case of coherent states which saturate the uncertainty relations. 
For the harmonic oscillator, such states are squeezed Gaussian states of the form ip(q) = 
exp(—Ziq 2 + z 2 q + z 3 ) with three complex numbers Z{ such that Re^i > 0. While Rez 3 is 
fixed by normalization and Imz 3 is only a phase factor, z\ — a.\ + if3\ and ^2 = ^2 + iftz 
determine the peak and fluctuations of the state. One can easily see that 

(?) = (46) 

<P> = 01102 ~ 01201 h, (47) 

G " = 2 (48) 

Gpp = ai h 2 + ^h 2 , (49) 



n 

2ori 



G qp = -p-h (50) 



such that indeed 



G qq G pp - {G qp f = hi 2 . (51) 

As is well known, these states even describe dynamical coherent states of the harmonic 
oscillator, i.e. their form is preserved during evolution. Saturation of the uncertainty 
relation leaves two free second moments such as one spread parameter G qq and squeezing 
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G qp . Unsqueezed states imply (5\ — = G qp and only one free parameter specifies the 
width of the Gaussian. 

Our system is different and no operator for c (which would be an analog of q) exists. We 
have to work with exponentials instead, and thus even kinematical coherent states change 
compared to the harmonic oscillator. For dynamical coherent states the form must anyway 
be different since we are not dealing with the harmonic oscillator Hamiltonian. Thanks to 
the solvability of our model we are still able to determine properties of dynamical coherent 
states explicitly. This provides an instructive example of how properties of coherent states, 
and physical implications, can change when a system is not closely related to a harmonic 
oscillator. 

We thus look at saturation of our uncertainty relations (1431) . (jUJ) and (145]) . first remov- 
ing the 0-dependence. Subtracting (l4"5l) divided by four from (JHJ) yields 

c 2 - c\ - 4c 5 c 6 = -h 2 AB . (52) 

The coefficients of e~ 2 ^ and e 2 * in (jHJ) then yield two more independent relations 

(ci-c 2 )c 3 -4 = h 2 A\ (53) 
(ci - c 2 )c 4 - c 2 = hi 2 B 2 (54) 
in addition to (1431) which becomes 

4c 3 c 4 - -(ci + c 2 ) 2 = h 2 H 2 . (55) 

There are thus four equations for six variables, such that two remain free as in the case 
of Gaussians. Note, however, that the reality condition relates one of these, c±, to the 
expectation values, or the constants A, B and H . Thus, only one combination of the 
uncertainty parameters is free, which we can take as C\ — c 2 = 2(AH) 2 , and saturated 
states are more restricted than for the harmonic oscillator. (This is a consequence of the 
non-linear reality condition which relates some fluctuations to expectation values.) 

Without loss of generality we can assume A = B for bounces or A = —B for non- 
bouncing solutions since it simply amounts to choosing the origin of time such that the 
bounce (or the transition through p = 0) occurs at <fi = 0. Subtracting fl53|) and fl54|) then 
gives 

(ci - c 2 )(c 3 - c 4 ) = Cg - c 2 (56) 

which shows that ^-fluctuations are the same before and after the bounce, i.e. c 3 = c 4 , if 
and only if also [ cq | = \cq\ (recall that C\ — c 2 = 2(AH) 2 cannot be zero). This case will be 
discussed below. 

Using A = ±B, we can combine fl52l . fl53l and (|5~4"1) to obtain a further ^-independent 
relation 

(ci - c 2 )(c 3 + c 4 ) - (cs ± c 6 ) 2 ± i(c 2 - c{) = . (57) 
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Solving for C3 + c 4 and combining it with (1561) gives 

c 3 = ±Cs + C6 Cg j- l( Cl + C2 ) (58) 
ci - c 2 4 
c 5 ±c 6 1 

c 4 = c 5 ± -(ci + c 2 ) (59) 

ci - c 2 4 

in terms of only c\, c 2 , C5 and Cg. Using this in (1551) and combining it with (1521) . we have 
h 2 H 2 + i( Cl + C2 ) 2 = 4c 3 c 4 

/ C 5 ±Cg \ 2 Ci + C 2 2 1 2 

= ±4 c 5 c 6 ± (c 5 ±c 6 ) +-(ci + c 2 ) 

V Ci — C 2 / Ci — C 2 4 

= ± f C -^) 2 (4c 5 c 6 + c? - c 2 2 ) + h Cl + c 2 f 
\C1-C2J 4 

= h 2 A 2 ( C -^X + ](c 1 + c 2 ) 2 (60) 
and thus (with c x - c 2 = 2(AH) 2 > 0) 

|c 5 ±c 6 | = — (ci-02). (61) 
This shows that it is impossible to have both C5 and Cq zero, i.e. 

Q P ,i(j-j) ^ q and there 

arc 

always correlations between p and H which evolve in time. Together with (1521 . (|58|) and 
(!59j) . this last relation allows us immediately to express all parameters in terms of only c\ 
and c 2 , and all relations for saturation are solved. 

We now ask whether it is possible to have a coherent state which behaves semiclassically 
at one (late) time and has identical fluctuations before and after the bounce. We thus focus 
on the case C5 = ±c$. This assumption allows us to solve flBTl) directly for C5, 

M = |c 6 | = 2^(ci -c 2 ) (62) 



and to insert it in (|58|) . 



F 2 1 

cs = c 4 = ^42 (ci -c 2 ) ± ^(ci + c 2 ) . (63) 



Consistency with (Iq"5T) implies 

4^ = /l 2 ^ 2 + i(ci + c 2 ) 2 = (^( Cl - c 2 ) ± i(d + c 2 )) 2 . (64) 

Coherent states with identical fluctuations before and after the bounce are thus possible if 
and only if 

^(c 1 -c 2 ) 2 T^(cl-c 2 1 ) = h 2 H 2 . (65) 
17 



Using Cl - c 2 = 2{AH) 2 and Ci + c 2 = -2{AH) 2 + 2c x = -2(AH) 2 - 2H 2 ± 2A 2 + \fi 2 
(imposing the reality condition (I35j) ). we must then solve 

H 2 =p A 2 t A rr . 4 A 2 (# 2 =F A 2 — 1ft 2 ) A 4 ft 2 



, (Affh- — ^-—(AH 2 )-^- — 

H 2 y ' ^ H 2 y 1 H 2 4 



giving 



A 2 (H 2 iA 2 - ±h 2 ) \ A\H 2 tA 2 - \h 2 ) 2 + A\H 2 t A 2 )h 2 
1 J 2{H 2 tA 2 ) 2{H 2 tA 2 ) 
T A 2 ({H 2 TA 2 -\h 2 ) + a\H 2 T A 2 + \h 2 \) 
= 2{H 2 T A 2 ) • (66) 

Here we distinguished the two roots of the quadratic equation fl65|) by o = ±1 since another 
± has already been used for the two cases A = ±B. 

Depending on the signs involved there are four possibilities to have positive solutions 
for (AH) 2 : 

1. A = B, in which case there is a further distinction 

(a) A 2 < H 2 + \h 2 : Only o = 1 is allowed, implying (AH) 2 = A 2 . 



(b) A 2 > H 2 + jh 2 : Both signs for a are allowed, 

i. a = — 1 implies (AH) 2 = A 2 as above; 

ii. cr = 1 implies 



A 2 h 2 

(^) 2 - < 67 > 



2. A = — B, which allows only one choice of signs for a positive 



A 2 h 2 



There are two cases where (AH) 2 = A 2 which can satisfy the basic condition (AH) 2 H 2 
only if A 2 H 2 . Thus ci must be of the order H 2 . This can only happen if the bounce 
scale p(0) = A — |ft is small compared to the total energy, i.e. the universe enters the deep 
Planck regime during the bounce. However, large c\ (and thus large c 2 since C\ — c 2 must 
remain small) of the order H implies, using ( 163ft . that C3 is of the order H, too, and Ap 
is not small compared to p (in fact, not even smaller). This case does not give rise to 
semiclassical states at any time. 

The case A = —B allows small AH. However, as we already saw, the reality condition 
allows such solutions only if c\ = —H 2 —A 2 +\h 2 is large and negative. For states saturating 
the uncertainty relations, this implies that C3 is of the same size as —c\ and thus too large 
for the state to be semiclassical at any time. 
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s 

Figure 2: Two bouncing solutions for the expectation value of p and the spread around 
it. Generic states have different spread before and after the bounce (dashed lines), while 
unsqueezed Gaussian initial states lead to solutions which are symmetric around the bounce 
not only in their expectation values but also in spreads (solid lines). 



For the last possibility (1671) . (AH) 2 can be small compared to H 2 . For instance, if C\ is 
of the order Ah, which is allowed for semiclassical states, we have (AH)/H ~ V 'Ah/H ~ 
a/ h/H. Moreover, (Ap)/p ~ \Jh/A such that this final possibility does allow semiclassical 
states with equal spread before and after the bounce. This confirms the earlier indication 
that equally distributed fluctuations typically satisfy C3 ~ HH. Although there is a factor 
of h, fluctuations are rather large due to the factor of H which, for a universe with large 
matter energy, is a large number. It is possible to have smaller fluctuations which are not 
magnified by the matter energy, but only at one side of the bounce and at the expense of 
having much larger fluctuations at the other side of the bounce. 

There is an easier way to have symmetric fluctuations if one does not require that all 
uncertainty relations be saturated. One can argue that (jUj) is of primary interest since it 
determines the fluctuations of p and H, and that only this relation should be saturated. 
If this is done, symmetric fluctuations before and after the bounce are easily allowed. 
However, this can only be put in by assumption and not be inferred from conditions at 
one time after the bounce: only one of the relevant parameters C3 or C4 is controlled by the 
uncertainty relation at one late or early time <fi while the other one would be suppressed 
exponentially by e -2 '^'. Symmetric fluctuations before and after the bounce can thus not 
be proven but only be assumed for coherent states of this system. Generically, even a 
universe which is semiclassical at late times can have been highly quantum before the big 
bang. Examples for one symmetric and the generic non-symmetric bounces are shown in 

Fig. m 
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3.4.2 Gaussian states 



This seems to be in conflict with recent numerical results (35J [15J [32] where the p-fluctuations 
turned out to be very close to each other before and after the bounce. Indeed, this gave 
rise to statements to the extent that the universe was as semiclassical before the big bang 
as it is now. To resolve the apparent contradiction we have to look at unsqueezed Gaus- 
sian states as they were used by construction in those numerical simulations (explicitly 
removing phase factors so as to de-squeeze states), but are not covered by our preceding 
coherent state analysis. As emphasized before Gaussian states do not saturate the uncer- 
tainty relations in this system. They are thus not coherent, and not distinguished as they 
would be for usual quantum mechanical systems. Nevertheless, their properties are quite 
interesting. 

Let us thus assume that we have a state of the form ip(p) = iVexp(— zip 2 + z 2 p), 
supported on integer p, with Kezi > 0. Although these states have the same form as 
in standard quantum mechanics, the representation of basic operators is different. For 
instance, we have 

(P> = 5>l^)! 2 ' ( 69 ) 
p 

(j) = ^^(p)^(p-n) (70) 
p 

where we sum over integers and the shift by h arises from the action of e lc which is a 
shift operator in p. For states which are nearly constant on the discrete scale of p, the 
expressions can be seen as Riemann sums and approximated by Gaussian integrals which 
one can compute explicitly. 
We thus obtain 

P « ^, (71) 

J w - + + l(3lH exp(-(a 2 + A 2 )ft 2 /2tti ~ i(Jh ~ e*2/y«i)ft) , (72) 

G pp « (73) 
4a:i 



G pJ « fU 2 + JpL [h + -M exp(-(a 2 + P 2 1 )h 2 /2a 1 - ih{(3 2 - a 2 /V«i)) (74) 



1 01 2 ipxh ( a 2 

A + ~A 

4«! Aa{ 4«i \ a.\ 



G JJ ^(exp(-(a 2 + /5 1 2 )/l 2 /al)( 7 ^ + ^-^ 2 + ^ + ^^l (To) 



4cti ' Aa'f d{ ~ ' 2«i ' a'f J Aa'f 4 
1 ;h 2 — — i^—{h + ati/oci) J exp(— (a 2 + f3 2 )h 2 fat\ — 2i(f3 2 — a 2 f3i/ai)h) 



Aa{ 2ai 2a>i 



« G JJ 

„jj 1 ft 2 a 2 h a\ ( a\ a 2 h h 2 (3?h 2 . ., .., , ,, 

K 4^ + T + 2^ + 4"(4 + ^ + T + %l exp( - (aJ + ft),i/a ' ) - (76) 
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These expressions involving J are much more messy than the corresponding ones for Gaus- 
sian states in standard quantum mechanics. This is a consequence of the fact that Gaussian 
states for the system considered here are not natural at all. 

Nevertheless, the expressions simplify somewhat if one assumes that the state is un- 
squeezed, Pi = 0. Moreover, for semiclassical states we can use p 3> Ap 3> h~ which implies 
a 2 /ai ^> a 1 ^> h. Then, the leading order contribution to the energy fluctuations, 
derived using 

Re G JJ = ((^ + ^ + ^n)e-^ 2 -^ 2 -^-^h)e~^ 2 cos(2p 2 h) 



4a i 4af 2ai J 4af 4 2ai 



4a 



■h 2 cos(2p 2 h) 



and 



4ai 2 2ai 4af \4af 2a x 4 



4ai 



is 



G HH = --(G JJ + G JJ ) + \g jJ « ^-^ 2 cos 2 (/? 2 n) (77) 



4 V '2 4ai 

and with if = ImJ « — a 2 sm{p 2 K)/2ai we have 

QHH 



H 2 



axh l cot\p 2 K) (78) 



whose left hand side is constant throughout the whole evolution. At late and early times 
P 2 becomes small such that ai oc cot -2 \P 2 K) ~ sin 2 (/3 2 ^). This relates ai which determines 
the p-fluctuation to p 2 which determines the peak position of the wave packet. Thus, for 
unsqueezed Gaussians we prove that the fluctuations before and after the bounce are the 
same at times where sin 2 (/32^) takes the same value. Since the expectation value solutions 
(fl3|) are symmetric around the bounce point for any state, fluctuations of unsqueezed 
Gaussian states are shown to be symmetric around the bounce. 

This reconciles our calculations with the numerical calculations of [13, [32] and reinforces 
their validity. However, it also demonstrates that the result of identical fluctuations before 
and after the bounce is a consequence not of the generic dynamics of semiclassical states, 
but relies on the assumption that states are unsqueezed Gaussians. It is then not very 
surprising to find symmetric spreads since there is a single parameter determining the 
state, other than its expectation values. The fact that G HH is constant then implies 
directly that there is a fixed relation between this parameter, a±, and the peak position. 
Since expectation values are symmetric around the bounce for any solution, the spread 
must satisfy the same symmetry in this restricted case. As discussed before, there is no 
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intrinsic basis in this model to restrict states to such a form. They do not saturate the 
uncertainty relations, and even if the Gaussian form is assumed but general squeezing is 
allowed fluctuations before and after the bounce become independent of each other. There 
is then an additional parameter ft and only a certain function of spread a± and squeezing 
ft is fixed at the bounce-reflected point. This does not suffice to fix the spread to be 
symmetric. The precise relation follows by estimating G HH as above, now keeping ft / 0. 
Then, G HH /H 2 = const with 

Q HH ^ "I (q 2 + p2 )h 2 c()s2((/52 _ a2f3l/ai)h) 

and 

OL 2 ft 

H ^ — - — sin((ft - a 2 p 1 /a 1 )h) + - — ftcos((ft - a 2 Pi/a 1 )h) 
lct\ Lct\ 

provides the relation between a,\ and ft in terms of the expectation values. 



3.4.3 The role of superposed branches 

One could suspect that non-symmetric fluctuations are a consequence of small admixtures 
of superposed expanding branches on an initial contracting one, while superpositions of 
energy eigenstates of definite sign might always have symmetric fluctuations. As illustrated 
in Fig. El we choose two states and \ip 2 ) at (f> = — </> which are peaked at the same value 
p of p with the same p-fluctuations. If the first state is expanding, (tpi\H\tpi) =: H > 0, 
while the second one is contracting, (i/j 2 \H\ip2) — ~('4 J i\H\ipi) < 0, they are sharply peaked 
at different if, and thus (V'llV^) ^ 1- Any superposition \vp) = (\ipi) + a\tp 2 ))/ yl + \ a \ 2 , 
where \ifj 2 ) presents an admixture of a negative- if state, is then sharply peaked at Pq. At 
<p = <fto, expectation values of p in the two states of the superpositon have evolved away 
from each other and the combined spread Ap can be much larger than individual spreads 
A\p and A 2 p measured in \ipi) and \i/j 2 ), respectively. 

This can be analyzed more quantitatively: For simplicity, we assume that the H- 
fluctuations of the two states are nearly equal, AiH w A 2 H. In the state \ip), we have 



(AH) 2 = ^|^ 2 |^)-^|^) 2 ^ T -V^((^ 2 )i + l«| 2 (^ 2 ) 2 )-^ 2 f^ 

1 + \a\ \1 + \a 

' r((Aii/) 2 (l + M 2 ) + (A 2 H) 2 \a\ 2 (l + \a\ 2 ) + A\a\ 2 H 2 ) 



2\ 2 



2\2 



where subscripts of 1 and 2 at right brackets and A indicate which state is used for the 
expectation values. Now assuming (AiH) 2 w (A 2 H) 2 , we obtain 



(AH) 2 ^ (A 1 Hf 



H 2 H 2 



1 A 


a 2 \ 




"'(1 + 1 


a 


2 ) 2 



(79) 



Thus, the admixture changes the if -fluctuations only slightly and preserves Aii <C H if 
\a\ is sufficiently small. It would thus be allowed in our approximation, although not in a 



22 



1 1 1 ' 

-<t>o <t>o 

Figure 3: Sketch of two wave functions with the same peak position p at — (f>o, but one being 
collapsing the other expanding. After the bounce of the first solution at = 0, the wave 
packets deviate strongly at O - Both wave packets are illustrated by their expectation values 
and spreads, assumed to be symmetric around the bounce for the sake of the argument. 



superposition of only positive energy eigenstates. We will now see how the admixture can 
influence spreads of p before and after the bounce. 

While (AH)/H is constant in time, (Ap)/p changes and the initial value which by 
construction is close to (Aip)/p is not conserved. Using the behavior of exact solutions 
of expectation values, one can compute p-fluctuations at the bounce-reflected point of the 
initial state. At this time, the two states will, in the ^-representation, have evolved away 
from each other since one state corresponds to an expanding branch and the other to a 
collapsing one. We have 



(Ap) 2 = — V^((p 2 )i + |«| 2 (p 2 ) 2 )- 
1 + a 



(p)l + l 


a\ 


2 (Ph 


1 + 


a 





1 



\2\2 



(((A lP ) 2 + M 2 (A 2 p) 2 )(l + H 2 ) + M 2 0S> 2 + M 2 0S> 2 - 2|a| 2 (p) 1 (p) 2 ) 



(1 + \a\ 

= (AlP) + (1 + H 2 )2 ((P>1-(P>2) • 

In the last step we used the fact that by construction j^) is in the expanding branch at all 
times considered, thus (Ap)/p m const away from the bounce and (A 2 p) 2 ~ (p)l(Aip) 2 / (p)\ > 
(Aip) 2 . For instance, at the bounce-reflected O of j^i) with large O , we have 

(p> 2 (0o) ~Poe 2(t>0 » (p)i(0 o )(2po/#) 2 
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and thus 



(Ap)(0o) 2 = (A lP )(<l> ) 2 l + A\a\ 2 p 2 /H 2 16H 2 (Poy(j>hW (m 
p{<Po) 2 p(4>o) 2 1 + M 2 (l + |a| 2 ) 2 \h) p(0 o ) 2 ' 1 } 

Due to the large factor p\jH 2 ^> 1 for an initial state peaked at large volume, the p- 
ffuctuation have grown much more than the if-fluctuation if we pass through the bounce. 

To verify that this implies non-symmetric fluctuations we have to compute Ap at the 
bounce-reflected point of (p), which is not 0o due to the contribution from \ip2)- From 



(P) = 



(p)l + 1 


a 


2 (Ph 


1 + 


a 


2 



w (p)i + |a| 2 (p) 2 



for small |a|, we have 

(p}(0) « Fcosh0+ ^|a| 2 ife* +200 

which has its minimum at ^bounce = — log \Jl + |a| 2 e 2<?;>0 . The bounce-reflected point for 
(p) is thus at ^reflected = ^bounce + (^bounce - ( -<j>o) ) = 0o - log( 1 + | a | 2 e 2<t>0 ) . Evaluating (JHOJ) 
at this point, for simplicity assuming \a\ 2 e 2< ^° ^> 1 which presents a characteristic example 
for an admixture at large initial volume, gives 

(Ap) (0 reflected ) 2 _ (A P )(-0 O ) 2 1 | q | 4 + | Q | ^e~ 2 ^ 1 

^reflected) 2 ~ P(-<A)) 2 |«| 2 | « | 4 + 2e~ 2 ^ + H"^" 2 ^ |tt| 2 (l + |tt| 2 ) 2 ' 

in terms of the initial values at — 0o- The additional inverse powers of the small |a| com- 
pared to the spread of H show the growth of p-fluctuations after the bounce, demonstrating 
that the superposition will not have symmetric spread. 

One could thus suspect that such an admixture, which would not violate our condition 
AH <C H if a is sufficiently small, could be the reason for unequal dispersions before and 
after the bounce, while solutions of exactly positive H would have symmetric dispersions as 
explicitly shown for unsqueezed Gaussians. This conjecture cannot be true, however, be- 
cause the solutions we studied earlier only refer to expectation values and dispersions and, 
due to the decoupling in our solvable model, are completely independent of higher moments. 
The preceding construction of the admixture does provide states with suitable initial dis- 
persions and expectation values, but the specific states \ip) = (|V>i) + a\ip 2 ))/ \/l + \a\ 2 
also have fixed higher moments. There are many other states having the same expecta- 
tion values and dispersions but different higher moments, not corresponding to what one 
obtains from an admixture of a negative energy state. Such states allow for positivity as 
well as non-symmetric dispersions, which thus cannot be an artifact of a negative energy 
admixture. In short, the calculation confirms the intuitive expectation that an admixture 
does give non-symmetric spreads even if each state in the superposition has symmetric 
spreads. But it does not show the converse, namely that non-symmetric fluctuations could 
only be caused by an admixture of a negative energy state. 
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4 Conclusions 



A solvable model such as the one discussed here allows a detailed analysis of dynamical 
coherent states which would otherwise be difficult to handle. This provides valuable in- 
formation for quantum cosmology, just as the harmonic oscillator does traditionally for 
quantum optics (see, e.g., [36]). Compared to the harmonic oscillator, our system shows 
several new properties with implications for cosmology. 

4.1 Spreading of states 

Although the spreads are not constant for solutions to our system, ratios such as (Ap)/p 
are nearly constant in each pre- and post-bounce branch. Nevertheless, this ratio can, and 
in general does, take different values in both branches. Dynamical coherent states which 
exactly saturate the uncertainty relations and have spreads symmetric around the bounce 
point do exist but are not generic. Nevertheless, they are distinguished in a certain sense 
and can thus be seen as analogs of the harmonic oscillator ground state (although no unique 
symmetric state exists). Indeed, for such states the scale of fluctuations, (Ap)/p ~ y/K/H, 
is determined more sharply than without the symmetry assumption. 

Dynamical coherent states for the loop quantization are not Gaussians which turn out 
to have different and rather special properties. Nevertheless, Gaussian states can also 
be analyzed straightforwardly in this setting, with results being ultimately in agreement 
with recent numerical investigations. Interestingly, unsqueezed Gaussian states do have 
identical spread before and after the bounce. They cannot be coherent but may well serve 
as a special version of semiclassical states. This illustrates how differently coherent states 
in a new system can behave from those well-known for the harmonic oscillator. It also 
shows that a dynamical analysis of coherent and semiclassical states is always necessary 
even to select suitable initial states to be evolved. While dynamical coherent states are 
difficult to describe in most systems, solvable models make this possible which is now also 
available for cosmological systems. 

As demonstrated, the techniques of [5] reviewed in Sec. [2] provide an efficient way 
to derive coherent state properties. Moreover, they allow one to see how properties can 
change under perturbations away from the solvable model, such as by including a scalar 
potential. Although this has not been dealt with in the present paper, our detailed analysis 
of dynamical coherent states of the solvable system provides crucial information for the 
zeroth order of such a perturbation theory. 

4.2 Applications to cosmology 

There is only scant information on properties of the present state of our universe other 
than that it is, to a high degree, semiclassical. But this does not tell us which semiclassical 
or even coherent state describes it best. From the harmonic oscillator or free quantum 
field theories one is used to unsqueezed Gaussians to represent the vacuum state. But the 
form of such states depends on the system being looked at, and quantum cosmology is 
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not close to either the harmonic oscillator or free quantum field theory. Moreover, there is 
no obvious vacuum state for quantum cosmology, and even if a Hamiltonian should allow 
a ground state, it is unlikely to describe a universe able to expand to large volume. A 
dedicated analysis as done here shows which coherent states are available and what their 
generic properties are. These properties must be taken into account for robust statements 
about quantum cosmological systems. 

For instance, we have seen that squeezing must be allowed for generic states, which 
crucially changes properties such as the symmetry of spreads before and after a bounce. 
In fact, squeezing of semiclassical states often plays a large role for decoherence or the 
transition to classical behavior (e.g. in the context of inflation [37J EH [39]). Thus, a 
semiclassical state at large volume of a cosmological model should indeed be assumed 
to be highly squeezed. Lacking additional input, robust cosmological conclusions can be 
drawn only with reference to generic coherent states. Then, no strong restrictions on 
fluctuations of a state before the bounce are justified. For all we know, it could have been 
coherent but with large quantum fluctuations. 

4.3 Effective equations and the possibility not to bounce 

Although we have proven that the solvable quantum system is exactly described by the 
effective Hamiltonian 

H = (H) = ~ J) =psinc 

(determining the same equations of motion for p and J as H does for (ft) and (J)) it 
is possible, depending on the initial state, that the system does not bounce but reaches 
p = in finite proper time. This looks contradictory at first sight since such an effective 
Hamiltonian implies the effective constraint equation 

tj2 2-2 2 

H = p sin c oc . 

With p^ constant and | sin c\ < 1 there must be a non-zero lower bound for p, the bounce 
scale. 

However, while the effective Hamiltonian of our linear system does not receive correc- 
tions from quantum back-reaction of fluctuations and other quantum variables, the reality 
condition JJ^ = p 2 is non-linear. Classically, this condition implies that c is real and thus 
| sine] < 1. But the reality condition fT2oT) now receives quantum corrections of second 
order in quantum variables which, for suitable states, can remove the bounce. A state can 
then enter the "classically forbidden region" where |sinc| > 1 while still respecting the 
quantum reality condition. 

As shown explicitly, zero volume is reached when the parameter C\ is negative and 
large. While this can be achieved respecting the spositivity condition AH <C H, such 
a state would never be semiclassical. Thus, any state which is semiclassical at one time 
will give rise to a bounce. These are all relevant states since a boundary condition for 
modelling our universe, however distantly, is always that there is at least one large volume 
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regime in which the state is semiclassical. But the possibility of states which do not 
bounce demonstrates the non-triviality of the result. Simply replacing d 2 in the Friedmann 
equation by a bounded function is not enough; any such replacement would have to be 
followed up by a coherent state analysis which is much more non-trivial than an analysis 
of the resulting "effective" Friedmann equation obtained by the naive replacement. Loop 
quantum cosmology with a free scalar passes this more stringent test and thus provides 
the first example in loop quantum gravity where complete effective equations have been 
computed. 

We end by repeating that any physical statements derived from a single model have 
to be confirmed by a perturbation analysis around the model. This is feasible in our 
case, as it is for perturbations around any solvable model, but still requires detailed work 
which is now in progress. Only such an analysis could justify the transfer of results from 
single models to our own universe. It may well be that this removes the bounce through 
back-reaction of quantum variables G a,n on the expectation values. In particular, it is 
then conceivable that a state starts out perfectly semiclassically at large volume where its 
expected volume collapses, evolves for a long time to small volume and all along picks up 
corrections from quantum back-reaction. Since also quantum variables evolve, it cannot 
be ruled out without further analysis that the analog of C\ does become negative and large 
close to the would-be bounce. If this happens, the bounce is avoided for the self-interacting 
state even if it starts out semiclassically by all possible conditions one could pose. This is 
only one possibility out of many which can only be ruled out by performing a comprehensive 
perturbation analysis for which the results of this paper present the zeroth order basis. 
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